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ABSTRACT 

We present the first results on the history of star formation in the Universe based on the 'cosmic 
spectrum', in particular, the volume-averaged, luminosity-weighted, stellar absorption line spectrum of 
present day galaxies from the 2dF Galaxy Redshift Survey. This method is novel in that unlike previous 
studies it is not an estimator based on total luminosity density. The cosmic spectrum is fitted with 
models of population synthesis, tracing the history of star formation prior to the epoch of the observed 
galaxies, using a method we have developed which decouples continuum and spectral-line variations and 
is robust against spectrophotometric uncertainties. The cosmic spectrum can only be fitted with models 
incorporating chemical evolution and indicates there was a peak of star-formation rate in the past of at 
least three times the current value and that the increase back to z = 1, assuming it scales as (1 + z)^ , 
has a strong upper limit of /3 < 5. We find in the general case there is some model degeneracy between 
star formation at low and high redshift. However, if we incorporate previous work on star formation at 
z < I we can put strong upper limits on the star- formation rate at z > 1: e.g., if /3 > 2 then the SFR for 
1 < z < 5 scales as (1 -I- z)" with a < 2. This is equivalent to stating that no more than 80% of stars in 
the Universe formed at z > 1. Our results are consistent with the best-fit results from compilations of 
cosmic SFR estimates based on UV luminosity density, which give 1.8 < j3 < 2.9 and —1.0 < a < 0.7, 
and are also consistent with estimates of figtars based on the K-hand luminosity density. 

Subject headings: cosmology: miscellaneous, observations - stars: formation 



1. INTRODUCTION 

The analysis of the comoving star-formation rate (SFR) 
density as a function of redshift has been the subject of 
much recent work. The onset of large redshift surveys at 
z < 1 (e.g. Lilly et al. 1996) and z > 3 (e.g. Steidel et al. 
1999) has allowed the volumetric emission of luminosity in 
different bands to be traced with redshift. In particular 
from these studies there is now good evidence for a rise in 
star-formation rate by a factor of about 8 between z = 
and z = 1 (Hogg 2002) and either a z > 1 decline (Madau 
et al. 1996) or plateau (Pettini et al. 1998). 

Most measurements of the SFR to date have been based 
on some type of luminosity density which is thought on 
theoretical and/or empirical grounds to trace the star- 
formation rate. This use of luminosity per unit volume re- 
flects an attempt to decouple the stellar history of the Uni- 



verse from its dynamical history. Other statistical mea- 
sures such as object counts versus luminosity or redshift 
are affected by the changes in the number of galactic ob- 
jects by merging as well as evolution of stellar populations. 
In contrast, the change in the light budget per unit volume 
is only affected by the stellar production in that volume 
regardless of the changes in the number of objects. Thus, 
for example, it makes sense to compare the total produc- 
tion of stars with the total metal abundance today (Cowie 
et al. 1988). 

Popular tracers of star-formation rate include: the UV 
1500-3000A continuum (Connolly et al. 1997; Madau, 
Pozzetti, & Dickinson 1998); the radio continuum 
(Mobasher et al. 1999); emission in the Ha and H/3 lines 
(Glazebrook et al. 1999), and; other line emission such as 
Lya (Kudritzki et al. 2000) or Oil (Colless et al. 1990; 
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Hogg ct al. 1998). The far-IR thermal dust emission has 
also been used to trace star-formation rate (Hughes et al. 
1998) although so far without the benefit of redshift in- 
formation. All these measures have in common the use of 
some kind of luminosity per unit volume whose change is 
proportional to the star- formation rate (with corrections). 
The debate over the z > 1 slope reflects the uncertainty in 
the dust correction to the UV continuum measurements, 
which are the easiest to measure at high redshift but most 
affected by dust. The other indicators are harder to mea- 
sure and are affected by small number statistics. 

In this paper, we present new constraints on the history 
of star formation based on the ensemble stellar populations 
of relatively nearby present day galaxies (z < 0.3). The 
concept is to use the average spectrum of nearby galaxies 
to constrain the earlier history of star formation leading 
up to that stellar population. The average spectrum con- 
tains absorption features for stars of all ages and probes 
look-back times of 0.2 -10 Gyr. With the advent of the 
large galaxy redshift surveys, such as the Two-degree Field 
Galaxy Redshift Survey (2dFGRS) and the Sloan Digital 
Sky Survey main galaxy sample (SDSS-mgs), it has be- 
come feasible to combine the spectra of 10^ 10''' galaxies 
to form very high signal-to-noise intermediate resolution 
spectra that represent the average emission of the uni- 
verse at various redshifts (0.03 ^ z < 0.25). Effectively, 
the surveys can be regarded as having a series of apertures 
on the cosmic background emission rather than apertures 
on individual galaxies. 

What is novel about this method is that it is not based 
on any luminosity output with time of the Universe, un- 
like all the other indicators discussed above. The method 
uses an integral over the star-formation history rather 
than attempting to track the derivative (the SFR) and 
uses the whole visible spectrum at intermediate resolu- 
tion. Perhaps the nearest approach to this in the past 
has been the work of Abraham et al. (1999) who used the 
color distribution of z ^ 0.5 galaxies to derive their star- 
formation histories which were then combined to form a 
cosmic star- formation history. Recently, Hopkins, Irwin, 
& Connolly (2001) estimated the global SFR density from 
star-formation history measurements of the Local Group. 
The problem here is that the Local Group SFR may not 
represent the cosmic mean and is subject to large cosmic 
variations such as recent "mini-bursts" of star formation 
in the Milky Way. 

An advantage of this kind of 'fossil cosmology' approach 
over the direct measurement is a reduced sensitivity to 
extinction. Young stars are born in dusty regions, this 
plagues the direct measurement approach. When they age 
they migrate out of such regions and contribute to the 
older stellar populations we observe. 

In this paper, we describe the application of this method 
to 166 000 spectra in the redshift range 0.03 0.25 from the 
2dFGRS (CoUess et al. 2001). The plan of this paper is 
as follows. In Section 2 we describe the 2dFGRS data 
and our methods for combining the spectra. In Section 3 
we describe our analytic models for cosmic star-formation 
scenarios. In Section 4 we describe our fitting procedure 
and the best-fit models. In Section 5 we discuss the im- 
pact of possible biases on our results from aperture efl[ects 
and luminosity selection. Finally in Section 6 we give our 



conclusions. 

2. THE 2dFGRS data 

The 2dF Galaxy Redshift Survey is a magnitude-limited 
spectroscopic survey (Maddox et al. 1998; CoUess et al. 
1999, 2001) using the Anglo- AustraHan Observatory's 2dF 
facility which is capable of observing up to 400 galaxies si- 
multaneously (Taylor & Gray 1990; Lewis, Glazebrook, & 
Taylor 1998; Lewis ct al. 2002). The magnitude limit of 
the survey is an extinction-corrected bj of 19.45 selected 
from the Automated Plate Measuring (APM) galaxy cat- 
alogue (Maddox, Efstathiou, & Sutherland 1990a, 1996; 
Maddox et al. 1990b). By the end of the 2dFGRS sur- 
vey (2002 January), up to 250000 unique galaxy redshifts 
are expected to have been measured. The survey covers 
approximately 2000 deg^ of sky distributed between the 
NGP and SGP in high-galactic latitude fields. A full de- 
scription of the survey geometry is given by Colless et al. 
(2001). 

In 2001 June, about 173 000 unique galaxy redshifts had 
been measured. This is the sample used in the analysis 
presented in this paper. The spectra are observed through 
a fixed 2".l-fiber aperture and the wavelength coverage 
varies only slightly from observing run to run. consistently 
covering the range 3700-7860 A with a 2.1-pixel resolution 
FWHM of 9.0A. All the spectra in the survey have been 
eycballed and assigned a quality Q from 1 to 5 (Colless 
et al. 2001): 1, no identifiable redshift; 2, a possible red- 
shift; 3, a 90% rehable redshift; 4, a 99% rehable redshift; 
5, a 99% reliable redshift with a high-quality spectrum. 
The survey is considered to consist of those galaxy spec- 
tra with Q > 3 (approximately 92% of galaxies observed). 
These galaxies have a median redshift of about 0.11 and 
the typical spectral signal-to-noise ratio at the survey limit 
is about 10 per pixel. 

We construct our 'cosmic spectra' in redshift slices z — > 
z + Az by applying an instrument response correction, 
de-redshifting to the rest frame and summing up all the 
Q > 3 spectra in the interval Az. The galaxies are scaled 
to match their measured 6j luminosity by comparison with 
an integration of the bj filter curve over the mcasiired spec- 
trum (in the observed frame). The scaling allows for the 
fact that the fibres sample only a fraction of a galaxy's 
light and for extinction and exposure-time variations be- 
tween observations. A maximum scaling is allowed (as a 
function of apparent magnitude) to avoid adding excessive 
noise from poor-quality spectra. This scaling limit is only 
applied to 5% of the data. Finally the cosmic spectrum 
is normalized to a mean of unity over a set wavelength 
range to re-iterate, our method of analysis uses spec- 
tral features not absolute luminosities other than for the 
weighting of the galaxy spectra. This spectrum represents 
the spectral emission per unit volume in the z ^ z + Az 
interval down to the limiting magnitude of the survey. As 
a result of the 6j« 19.45 Hmit of the 2dFGRS, the cor- 
responding limiting absolute luminosity will be higher at 
higher redshifts. This exclusion of lower luminosity galax- 
ies is a possible source of bias. This is quantified in Sec- 
tion 5. 

The instrument response correction consists of apply- 
ing the average 2dF spectral response as given by Lewis 
et al. (2002). Only a simple scaling is applied for light lost 
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outside the fixed angular size fiber aperture which again 
is another potential source of bias to be discussed in Sec- 
tion 5. 

Examples of cosmic spectra, at a series of redshifts, are 
shown in Figure 1. As might be expected from a broad 
star-formation history, the average emission from the Uni- 
verse looks remarkably like an Sb-Sbc galaxy spectrum 
(Kennicutt 1992). 

3. COSMIC STAR-FORMATION SCENARIOS 

First we must decide how best to parameterize cosmic 
star-formation histories. We use two different families of 
parameterization. 

The first is a physically motivated parameterization, us- 
ing an infall model to describe star formation as a function 
of time and timescales. 

The second is a purely empirical model, which describes 
star formation as a series of power laws with redshift. This 
facilitates direct comparisons with the literature. 

For the cosmology we consider two types, both con- 
strained to be flat in accord with the Cosmic Microwave 
Background data (de Bernardis et al. 2000): 

1. The emerging standard cosmology (which we de- 
note CI), which has cosmological parameters Co = 
{h,nmo,^Ao) (0.70,0.3,0.7)^ These parameters 
are the current best constraints from a variety of 
observations (Silk 1999). 

2. A longer-age cosmology (which wo denote C2) with 
Co = (0.55,0.2,0.8). These parameters have been 
adjusted on the longer age side to the margins of 
consistency with modern limits. 

For galaxies at 2; = 0.1 with ^form = 5, these cosmologies 
give ages of 11.0 and 15.7 Gyr, respectively. 

Given a star-formation history and a cosmology we 
can compute spectra using standard evolutionary synthe- 
sis codes. We used the PEGASE code (Fioc & Rocca- 
Volmerange 1997)^. Both parameterizations require an 
initial-mass function (IMF) and we choose the Salpeter 
(1955) power law slope with stellar mass in the range 0.1- 
120 Mq. We have also investigated the Kennicutt (1983) 
IMF over the same mass range. We also have choices to 
make on metallicity evolution and dust extinction. For 
the former, the enrichment of the inter-stcUar medium is 
determined using the calculations of Woosley & Weaver 
(1995) within the PEGASE code. For the dust extinction, 
the prescription for an inclination-averaged disk geometry 
is used. However, the choice of extinction, E{B — y) ~ 
0.2 ± 0.1, makes negligible difference to the results in this 
paper for two reasons: (a) we are only determining the 
relative SFR as a function of time or redshift, i.e., not 
comparing predicted luminosity with luminosity density, 
and; (b) the robust results are primarily constrained us- 
ing high-pass filtered spectra which arc insensitive to the 
extinction model. We note that there will be a second- 
order effect from any systematic variation of extinction 
with spectral type and/or luminosity. 

1 /i = Ho/100 km s-iMpc-i 

2 PEGASE Version 2, revised 2001 May 5, URL http://www.iap.fr/users/fioc/PEGASE.html 



3.1. Physical parameterization 

The first parameterization is a natural scenario. Star 
formation starts at z = Zform- Gas falls into sufficiently 
dense regions to form stars with an infall timescale ti'. 

Mgalaxy = 1 - C"*/*' (1) 

The gas in these regions forms stars at a rate (S) propor- 
tional to the amount of gas available with a star-formation 
timescale tg'. 

S = Mgas/ts (2) 

This parameterization is readily implemented in the 
PEGASE population synthesis code (Fioc & Rocca- 
Volmerange 1997). The code includes recycling of ejecta 
into the inter-stellar medium (ISM) and consistent evo- 
lution of the metallicity. An analytical approximation to 
the star-formation rate can be determined assuming an 
instantaneous-recycling approximation (/ is the mass frac- 
tion of stars that arc not returned to the ISM) . Using this 
prescription the star-formation rate evolves as 

For fti ^ tg, the solution for the star-formation rate with 
time is given by 



S = 



e-*/*i _ e-ft/ts 



(4) 



fti-tg 

with initial condition iS = at t = 0. The normalization is 
such that the total mass of gas available is unity. This pa- 
rameterization is useful in that it is physically motivated 
and that it allows for a rise at early times and a fall at 
late times of the universal star-formation rate. The de- 
viation between the approximation and the code becomes 
apparent when recycling from lower mass stars becomes 
significant. An example is shown in Figure 2 with / = 0.7 
and 0.8 for the recycling approximation. 

3.2. Empirical parameterization 

In the second parameterization we empirically make the 
star-formation rate a function of redshift. In particular, we 
make the star-formation rate between z = 5 (our choice 
of formation redshift) and z = 1 proportional to (1 -|- z)" 
and between z = 1 and z = proportional to (1 -|- z)'^, 
with matched SFR at z = 1. This form allows us to make 
a comparison of the best-fitting parameters with previous 
studies which measure star-formation rate with redshift 
via luminosity densities. In particular we can compare the 
value of f3 defined and measured by Hogg (2002) as well as 
other measurements of the variation of the comoving SFR 
density with redshift (e.g. Lilly et al. 1996; Madau et al. 
1996). 

The total amount of star formation (r) between z = 5 
and z = as a fraction of the mass available was normal- 
ized to a range of values from r = 0.3 to r = 1.4 (the SFR 
normalization). In other words, the mass of gas available 
for star formation is unity and the r is the total mass of 
stars formed since Zfomi- This can be greater than unity 
because of recycling of material back into the inter-stellar 
medium. Higher r values result in higher average metallic- 
ity because the fusion products released into the ISM by 
supernovae are more abundant relative to the remaining 
gas. 
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Fig. 1.— Averaged 2dFGRS spectra from various redshift bins (0.025-0.04, 0.06-0.07, 0.08-0.09, 0.10-0.11, 0.12-0.13, 0.155-0.175, 0.20- 
0.25). The spectra were normalized to a mean of 1.0 between 4200 and 5800A, and offset by 1.0 from each other. 
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Fig. 2. — Example star-formation scenario using an infall timescale and a star-formation timescale. The instantaneous recycling approx- 
imation, Equation 4, assumes a i? (= 1 — /) mass-fraction of stars formed are returned instantaneously to the inter-stellar medium. The 
normalized SFR is relative to the total mass of gas available for forming stars, the integral of the SFR can be greater than unity because of 
recycling of material back to the ISM. 



4. BEST-FIT STAR-FORMATION SCENARIOS 

4.1. Reductions 

The 2dFGRS data were divided into 14 redshift bins 
between z — 0.025 and z = 0.25, each containing about 
12000 spectra. Half the redshift bins were below z = 0.11 
(the median redshift of the survey). The lower and higher 
groups of 7 redshift bins were considered separately and 
together to test the robustness to varying aperture and 
selection effects (see Section 5 for a further discussion of 
this). 

For each redshift bin, the spectra were divided into ten 
positional bins based on their coordinates (4 regions in 
the NGP and 6 in the SGP). A normalized average spec- 
trum was calculated and the positional bins were used to 
estimate the uncertainties. This spectrum represents the 
total optical emission of all galaxies in the volume of the 
redshift shell down to the 2dFGRS magnitude limit. The 
absolute magnitude limit is fainter than or about M* out 
to a redshift of 0.2 (the depth is discussed further in Sec- 
tion 5). 

To flux-calibrate the spectra they were divided by 
the 2dF response function of Lewis et al. (2002) before 
coadding the fluxes at their rest wavelengths. Spectra con- 
tribute to the final spectrum in proportion to their 6j lumi- 
nosity. In addition, the averaged spectra were smoothed 
and resampled to match approximately the much lower 
20A resolution of the spectral library (Lejeune, Cuisinier, 
& Buser 1997) used by the population-synthesis code. 

4.2. Goodness of fit 

To evaluate the 'goodness of fit' between a star- 
formation scenario and the 2dFGRS data, we compare 
the spectrum from each redshift bin with the appropri- 
ate model spectrum (at the same age) from the scenario. 
Fiber spectra are known to be difficult to flux accurately 
so it is desirable to develop methods which are insensitive 
to small spectrophotometric uncertainties. 

The first method we used was to allow for the possibil- 
ity of spectrophotometric calibration errors by including 



a correction function before evaluating the fit. This was 
to account for spectrophotometric discrepancies between 
the 2dF response curve used to calculate the average spec- 
tra and the true average 2dFGRS response. We used a 
fourth order polynomial with observed wavelength for this 
spectrophotometric correction, the coefficients being de- 
termined by ratioing the model and data spectra (exclud- 
ing the strongest emission lines as per the fitting). An 
important point is that the correction is constrained to 
be the same function at all redshifts, thus the degeneracy 
between the spectrophotometric correction and the model 
fitting is partially broken because of the range of redshift 
covered (0.03-0.25). Typically in our fitting we find this 
spectrophotometric correction to be of the order 5-10% 
change in the response function of Lewis et al. for the 
best-fitting models (within 3a). These values represent 
the RMS relative difference over the wavelength range be- 
tween the polynomial correction and a constant value (rep- 
resenting changes in the relative, not absolute, spectropho- 
tometry). Since the response function was measured from 
standard stars observed during a single hour, such a dif- 
ference could arise from: (i) unaccounted for, reduction 
or throughput discrepancies between extended and point 
sources, and/or; (ii) variations in the instrument response 
over the 2dFGRS survey time. 

In evaluating the fit between a model spectrum and a 
2dFGRS data spectrum, we compare the low-pass and 
high-pass information separately to form two figures of 
merit (FOMs) from the normalized spectra. This is illus- 
trated in Figure 3. The normalization is to set a mean 
of unity over the rest-wavelength range 4100-6200A, with 
the lower limit being just above the 4000A break and the 
upper limit being just in the detected range for galaxies 
with z = 0.25. Both the normalized model spectrum and 
normalized data spectrum are smoothed using a top-hat 
function of width 200A (10 sampling points). The original 
spectra are divided by the smoothed spectra to produce 
the high-pass spectra. FOM A is the reduced value from 
the difference between the low-pass spectra and FOM B is 
is the reduced value from the high-pass spectra. 
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Fig. 3. — Example model and data speetra with the low-pass (A) and high-pass (B) information separated. The data have been adjusted 
for possible spectrophotometric errors as described in Section 4.2. FOM A is the reduced value determined from the difference between 
the smoothed spectra (offset by —0.3 in the plot). FOM B is determined from the high-pass spectra (offset by —1.0). Note that the strongest 
emission lines were excluded from the fitting procedure. The first plot shows a good fit with reduced values of order unity for both figures 
of merit. The second plot shows a poor fit for both FOMs (x^ jv ~ 17 and 6 for A and B, respectively). 



This high-pass FOM B is thus robust against any large- 
scale errors in the spectrophotometry and in fact we find 
FOM B is negligibly different whether we include a spec- 
trophotometric correction or not. FOM B places more 
reliable constraints on the best-fit star-formation scenar- 
ios since it is less affected by the systematic uncertainties 
of spectrophotometry and extinction. FOM A is an inde- 
pendent check: our aim is to get consistent histories from 
both the low-pass (continuum) and high-pass (absorption 
line) information. From the best-fit star formation we can 
obtain a range of probable spectrophotometric corrections 
that were used to determine the mean and uncertainty in 
the 6j k-corrections of the galaxy spectra (Madgwick et al. 
2002). 

In order to determine the uncertainties in our set of 
spectra we take the empirical approach of dividing up the 
survey into ten separate, approximately-equal sky areas. 
We compute the mean spectrum separately for each area 
and use the variance between them to set out errors as a 
function of wavelength. The RMS errors computed this 
way are around 0.2-1% for the high-pass spectrum and 1- 
3% for the low-pass spectrum. An additional uncertainty 
of 1.5% per wavelength bin is added in quadrature to ac- 
count for intrinsic model inaccuracies due to, for example, 
the signal-to-noise ratio of the spectra in the input atlas. 
This uncertainty value was chosen so that the reduced 
values for the best-fit star-formation scenarios were ap- 
proximately unity. 

For a given star-formation scenario, we average the 
FOM values over all the redshift bins. We calculate the 
figure-of-merit quantities by summing over all wavelengths 
except near strong nebular emission lines (Oil 3727; OIII 
5007; Ha 6563; Nil 6583; SII 6716 & 6730). We do this 
because our goal is to calculate the star-formation history 
from the stellar emission only. The nebular-line emission 
is a measure of the instantaneous star-formation rate in 
HII regions, but the PEGASE code uses over-simple as- 
sumptions to translate star-formation rate to a set of line 
ratios. Since understanding ionization levels and interpret- 



ing line ratios are complex subjects, and since the current 
strength of emission lines does not constrain past star for- 
mation rates, we choose not to use this information. In 
addition the emission lines will be contaminated by active 
galactic nuclei in a few percent of our galaxies, so ignoring 
the lines will minimize the bias this causes. 

4.3. Results 

For the fitting procedure we calculate a grid of model 
values, {ti,ts, Z{orm) for the physical parameterization and 
{a, P, r) for the empirical parameterization. We estimate 
that after our data processing there are approximately 
20 degrees of freedom for FOM A (about 22 independent 
wavelength points) and 200 degrees of freedom for FOM B. 
The minimum value was about unity for both figures 
of merit using the uncertainties described above. We be- 
lieve the additional uncertainty of 1.5% necessary to ob- 
tain a suitable value refiects the remaining imprecision 
of the model spectra, as these types of models were con- 
ceived originally to model galaxy colors. Improved, higher- 
resolution models will be required in future work. We 
proceed to interpret the star-formation laws of the best- 
fitting models. We note that the best fit matches the data 
to about 2% RMS, so in absolute terms we are accounting 
for the volumetric light emission very well. The problem 
remains that the high-pass data are much better than the 
current state of the art in population synthesis modeling. 

In order to estimate the confidence limits on the FOMs 
(Ax^/i^), we investigated the variation of the best-fit pa- 
rameters using Monte Carlo simulations. It is not sufficient 
to apply the standard confidence limits because of corre- 
lated uncertainties between different wavelength points. 
The investigation included using the spectra from the ten 
positional bins separately, applying random errors of 1.5%, 
using different redshift ranges (e.g., 0.03-0.11, 0.11-0.25) 
and varying the dust extinction in the models. As ex- 
pected, FOM B is fairly robust and we use Ax^/i^ = 0.15 
for the 3cr (99.73%) confidence limit. FOM A is strongly 
affected by using different positions on the sky and dif- 
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Fig. 4. — Best-fit regions of logts versus logt; for cosmology CI with Zform = 5. The two plots use different redshift ranges from the data. 
The contours represent the 2a and 3cr formal confidence limits, with dotted lines for FOM A and solid lines for FOM B. The FOMs were 
calculated using the average from 7 redshift bins at each grid point. The diamonds represent the FOM B best-fit parameters for each bin. 
Section 3.1 describes the cosmic SFR parameterization using a star-formation timescale (is) and an infall timescale (<;). 



ferent redshift ranges. This is principally due to spec- 
trophotometric uncertainties, and we use A^^/i/ = 2.1 for 
the 3(7 limit (when plotting contours using the full redshift 
range). In neither case is the FOM significantly affected by 
changes in the dust extinction used in the model. This is 
not surprising for FOM B since it is determined from high- 
pass spectra. A spectrophotometric adjustment is applied 
to the data before determining the FOMs and therefore 
FOM A is principally sensitive to variations over ranges 
of 200A (the smoothing length), such as the 4000A break, 
and is less sensitive to extinction. 

The best-fit parameters are fairly degenerate. Figure 4 
shows the best-fit regions in log ts versus log space with 
Ziarm = 5. The rcsults for the lower and higher redshift 
ranges are shown separately. The solid contours represent 
the 95.4% and 99.73% ('2cr' and '3cr') confidence bound- 
aries for FOM B, while the dotted contours are for FOM A. 
Notably, the regions contained by the FOM B boundaries 
at low and high redshift are similar whereas the FOM A 
boundaries are significantly different. Wc do not over- 
interpret this, noting that FOM A is less reliable due 
to spectrophotometric uncertainties. The best-fit model 
(FOM B) has U ^ 4000 Myr with U < 200 Myr (note de- 
generacy for low ti). Both redshift ranges have ti < tg 
within the 2a levels even though the SFR is similar on in- 
terchange of fti and ts (see approximation, Eqn. 4) . How- 
ever, the metallicity evolution is different depending on 
whether infall of new gas or the star-formation timescale 
determines the SFR at late times. 

The figure-of-merit values become degenerate for ti <C 
is (galaxies form quickly in comparison with the star- 
formation timescale) while maintaining a good fit. In Fig- 
ure 5, ti is set at 100 Myr and the best-fit regions of logtg 
versus logzform are identified for cosmologies CI and C2. 
At the 3(7 limits (FOM A and B) , the redshift of formation 



is greater than or about 0.65. 

To aid understanding of the meaning of our results. Fig- 
ure 6 shows a variety of star-formation histories that are 
within the 3cr limits. The first plot shows scenarios with 
ti ^ ts, i.e., identified from Figure 5 (CI). The second plot 
(Fig. 6) shows 'smoother' scenarios with ti ts, identified 
from plots similar to those shown in Figure 4, the right 
edges of the Sc boundaries. 

The best-fit scenarios cover a large range of Zform but 
have the following in common: the normalized SFR at 
z = is in the range 0.02-0.04 Gyr"! and 84-92% of stars 
formed prior to z = 0.3 for CI. 

It can be seen that the main constraint is on the slope 
of the star- formation rate for < z < 1. Fitting a SFR 
of (1 -I- z)'^ from z = to 1 to the natural scenarios 
with Zform > 1, we obtain from the 3(t confidence lim- 
its: 1 < /? < 4.5 (CI) and 1 < /3 < 4 (C2). If we fit a SFR 
of (1 -|- z)" from z = 1 to 4, there is no lower limit on a 
since at least one of the well-fitting models has little or no 
star formation prior to z = 1 . The upper limit from these 
natural scenarios is a < 1 for both cosmologies. 

The ii-tg-^form parameterization is limited in its scope 
for changes of SFR with time. Therefore, to further test 
star-formation history, we look at the a-^-r parameteri- 
zation. Figure 7 shows best-fit regions in a versus (3 with 
r = 1.1 for both cosmologies. 

In these scenarios, the galaxies start fully constituted 
with gas (there is no infall) and consistent evolution of the 
metallicity is implemented. Constant metallicity scenarios 
were also tested but found not to be consistent within the 
3(7 limit on FOM B. If we take a SFR normalization of 1.1, 
meaning a total mass of stars formed between z = 5 and 
z = equal to 1.1 times the mass of gas available, there 
is a degeneracy across the plane of a versus (3 (Figure 7) . 
Scenarios with /3 < cannot be ruled out for a > 2.5. 
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Fig. 6. — Example set of star-formation scenarios that are within Scr of the minimum FOM values for cosmology CI. The first plot shows 
scenarios with t; <C and the second plot with i; ~ tg. 



However, this would imply a minimum in the SFR around 
z = 1 which is in disagreement with many cosmic SFR 
density studies based on photometry. These also repre- 
sent a branch of solutions which are impossible to produce 
using the physical parameterization. The conservative up- 
per limit of this study is /3 < 5. Taking the conservative 
lower limit on /3 from the analysis by Hogg (2002), /3 > 1.3, 
we obtain a < 3 (for CI). If we assume a < (i.e., star 
formation declined or remained constant for z > I which 
corresponds to the lower-right branch of solutions in Fig- 
ure 7) then we obtain a range of 1.5 < /3 < 5. This is in 
very good agreement with Hogg's results which represent 
an ensemble of different luminosity-density based indica- 
tors. 

Figure 8 shows best-fit regions in a versus r with f3 — 3 
for CI. The two plots show the results for different IMFs. 
There is no lower limit on a given /3 = 3. There are up- 



per limits of a < 0.5 for the Salpeter IMF and a < 1 for 
Kennicutt IMF, with a < at the 2(T-confidence level for 
both IMFs. There is some degeneracy between the final 
metallicity {z — 0) and the IMF. For example, if we take 
star- formation scenarios within a > — 1 and the 3cr con- 
fidence boundaries, there are scenarios with around solar 
to twice-solar metallicity using the Salpeter IMF and with 
around half-solar to solar metallicity using the Kennicutt 
IMF. 

With /3 > 3, the best-fit models have a decrease in SFR 
prior to z = 1 {2a limit of FOM B). Figure 9 shows best-fit 
regions in a versus r with /3 = 2 for CI. Here, the best-fit 
models are consistent with a plateau in SFR prior to 2 = 1 
or a marginal increase or decrease. An upper limit would 
be a < 1.5 {2a), if /3 > 2 as is suggested by most direct 
methods of tracing the cosmic SFR (Hogg 2002). 
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Fig. 7. — Best-fit regions of a versus /3 for cosmologies CI and C2. See Figure 4 for contour meanings. Section 3.2 describes the cosmic 
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0.04 in the first plot, and; 0.01 and 0.02 in the second plot. 



5. POSSIBLE SELECTION EFFECTS AND BIASES 

Ah spectroscopic information is subject to the effect of 
aperture bias (e.g., see the discussion of Kochanek, Pahre, 
& Falco 2002). In particular, the aperture may be too 
small to encompass a representative fraction of a galaxies 
light. This effect will obviously be greater at lower red- 
shift. 

Figure 10 shows the increase in the effective size of a 
2dF 2".l-fiber aperture with redshift. At z = 0.1 this is 
2.7/i^^kpc which is comparable to a typical large disk 



scale length of 3 /i^^ kpc (de Jong & Lacey 2000). In prac- 
tice, the effective aperture is more like 3.4 kpc (2".6) 
as the median observation seeing was about 1".5. Thus it 
is reasonable that the effect of aperture bias will be much 
smaller for z > 0.1 as we are sampling more than half the 
total light of galaxies. 

The data allow for a further test of this because of its 
own internal seeing variations. Spectra taken in bad atmo- 
spheric seeing are a rough proxy for spectra taken through 
a larger aperture, as the object is smeared out over a disk 




Fig. 9. — Best-fit regions of ct versus r witli (3 = 2. 
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Fig. 10. — Selection and aperture effects in the 2dFGRS and, for comparison, the SDSS main galaxy sample. The left plot shows the 
limiting magnitude relative to the characteristic magnitude of the Schechter (1976) luminosity function versus redshift. For the 2dFGRS 
miim 19.45; m* -19.79 + SlogDt + 1.9z + 2.72^ (Madgwick et al. 2002); the solid lines show the limits with ±0.04 error in M* and 
±10% error in the k-correction. For the SDSS (r'): 17.7; m* ^ -20.83 + Slogi^L + 2.51og{l + z) (Blanton ct al. 2001); the dashed 

lines show the limits with ±0.03 error in M* and ±20% error in the k-correction. The luminosity distance is calculated assuming a world 
model with Co = (1.0,0.3,0.7). The right plot shows the aperture scale for fiber diameters of 2".l and 3".0. 



about the size of the seeing disk. As a test of this, we 
chose a sample of about 1500 galaxies (z < 0.15) that 
had measured spectra taken both in relatively good seeing 
1".5) and in poor seeing (> 3") such that the differ- 
ence in representative aperture was greater than a linear 
factor of 1.5 for each galaxy with Q > 3 for both spectra. 
With this sample, we measured the change in equivalent 
width (EW) of a number of lines between the larger rep- 
resentative aperture and the smaller aperture. In some 
individual spectra, the sky subtraction was inadequate for 
accurate EW measurement. Therefore, measurements for 
which the reduced counts were too low or contaminated 
by sky-emission lines were excluded. For further robust- 
ness, the galaxies were divided into subsamples of ten as 
a function of redshift and the median change in EW of 
each group was determined. The results are shown in Fig- 



ure 11. Similar results were obtained if a weighted mean 
(by counts) was used rather than a median. 

There is no evidence for any aperture effect in terms of 
the average change of EW for the galaxies in this sample 
given a difference in representative aperture with a linear 
factor in the range 1.5-2.3. We cannot rule out a small 
aperture effect on the measured cosmic spectra particu- 
larly for z < 0.05 where the test sample is small. 

How do the measured cosmic spectra vary as a function 
of redshift due to selection effects, aperture effects, if any, 
and evolution? Figure 12 plots the equivalent width of 
on, US, CH, OIII, Mgl and NaT as a function of redshift 
for a complete sample and for a volume-limited sample. 
These measurements were made on luminosity-weighted, 
averaged spectra determined from redshift bins of mostly 
0.005 in extent. To interpret this figure, we first consider 
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Fig. 11. — Change of line equivalent widths - On, Can, Mgi, Ha - between measurements taken in poor seeing and in good seeing. Each 
diamond represents the median change of a subsample of ten galaxies. The lines represent the mean of the subsamples, ±3x the standard-error 
of the mean, divided into three redshift ranges. 



the vohimc-limited sample (diamonds). 

1. In the redshift range 0.05-0.2, we expect aperture 
effects to be minimal. Therefore, the change in EW 
in this range could be primarily due to evolution 
and the changes agree qualitatively with a /? ~ 3 
scenario. It is not possible to make an accurate 
comparison due to the lower resolution of the model 
spectra. However, qualitatively they agree well: de- 
crease in nebular- line emission between z = 0.2 and 
0.05; increase in CH absorption EW; less increase in 
Mgl and Nal absorption and, not surprisingly, the 
H(5 absorption decreases due to a lower fractional 
contribution of A stars. The solid lines in Figure 12 
are fits to these changes: with the EW oc (1 + 2)^"^ 
for the emission lines, and EW = P'^^z + const, for 
the others. These constitute a detection of cosmic 
evolution. However, we cannot unambiguously sep- 
arate cosmic SFR evolution from galaxy-population 
evolution and so we have not used these measure- 
ments to constrain the cosmic SFR in this paper. 

2. In the redshift range 0.02-0.05, the change in EW 
for most of the lines levels off or reverses direction. 
As described above, we expect the evolutionary 
changes to continue from 2 = 0.1 to 0.02. Our inter- 
pretation is is that at low redshift the fiber aperture 
will only see a small part of the galaxy which for a 
blue-targeted sample is biased towards more recent 
star formation - thus, the nebular-line emission is 
higher and the CH absorption is lower. We note 
that for an individual extended spiral galaxy, a 2dF- 
fibre spectrum could be biased towards more quies- 
cent regions (e.g., the bulge) or more star- forming 
regions. Here, we are strictly refering to the average 
effect on the measured cosmic spectrum (the me- 
dian bias may still be towards the bulge). We note 
also that the nebular-line emission may be more 
sensitive to aperture and selection effects (and we 



ignore them in the fitting) than features at other 
wavelengths derived from older stellar populations. 

Considering the all-selection sample (crosses) from Fig- 
ure 12, it is clear that luminosity-selection bias is having a 
more significant effect than aperture bias, i.e., the differ- 
ence between not including and including the lower lumi- 
nosity galaxies in the cosmic spectrum is larger than the 
difference between using small and large physical apertures 
on a volume-limited sample given the above interpretation. 
This is most evident for the Oil and OlII emission lines 
showing that the less-luminous galaxies contain a higher 
fraction of recent star formation. The 2dFGRS is a magni- 
tude limited sample so the luminosity cutoff will rise with 
redshift. This is illustrated in Figure 10 - at = 0.1 we are 
probing 1.5 magnitudes fainter than the Schechter lumi- 
nosity function break M* . For a Schechter function with 
a faint end slope of —1.19 (the 2dFGRS value found by 
Madgwick et al. 2002) this will encompass about 70% of 
the light at 2 = 0.1, 95% at z = 0.04 but only 20% at 
z = 0.2. 

It should be noted that all magnitude-limited samples 
are dominated by M ~ M* galaxies. This is a conspiracy 
between the volume probed and the shape of the Schechter 
function (e.g., see Glazebrook 1992, for a discussion of 
this). Other studies of luminosity density evolution via 
flux-limited samples will also be affected by this domi- 
nance of M* galaxies, so the comparison of our results in 
Section 4.3 is consistent. 

If we consider our divided redshift ranges 0.03-0.11 and 
0.11- 0.25: for the lower range, the luminosity selections 
are less significant (the measured cosmic spectrum samples 
a significant range of the galaxy luminosity function) but 
the aperture effects may be biasing the averaged spectra 
towards more recent star formation; for the higher range, 
the aperture effects are less significant but the selection 
function is clearly missing a higher fraction of recent star 
formation from the lower luminosity galaxies. However, 
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Fig. 12. — Variation of average line equivalent widths - On, H(5, CH, Oiii, Mgi, Nai - as a function of redshift in the 2dFGRS data. The 
crosses represent 49 redshift bins between 0.020 and 0.32, selecting all galaxies. The dotted lines represent quadratic fits. The diamonds 
represent a volume-limited sample of 35 redshift bins between 0.025 and 0.20, selecting galaxies brighter than a k-corrected absolute magnitude 
of —21 (using h = 0.7). The solid lines represent straight-line fits over the redshift range where cosmic evolution is expected to dominate 
the measured variations. The normalized deviation is defined as (E — {E))/{E) where E is the equivalent width, and (E) is the mean of the 
complete sample. Note that the normalization is such that a positive value means an increase in emission for the emission lines (left panels) 
or an increase in absorption for the others. 



in general, we have found the results to be robust against 
choice of redshift range chosen: lower, higher or combined. 
In any case, the contours in Figures 4-9 take account of 
the variation in best-fit parameters between the lower and 
higher redshift ranges. 

6. DISCUSSION AND CONCLUSIONS 

We have developed a method of determining relative cos- 
mic star-formation history based only upon the spectral 
information. We find consistent results between both low- 
pass and high-pass spectral information and between phys- 
ical and empirical parameterizations of the star-formation 
law. If we assume the following: (a) the averaged 2dFGRS 
spectra (z in the range 0.03-0.25) represent the galaxy 
population as a whole; (b) the stellar populations can be 
approximated by one chemical evolution scenario, and; (c) 
the IMF, the PEGASE models and the cosmology are suf- 
ficiently accurate; then we can reach the following conclu- 
sions. 



1. 



The present day 'cosmic spectrum''^ is well deter- 
mined at z 0.1 and can only be fit with models 
incorporating consistent chemical evolution. Con- 



stant metallicity is strongly ruled out. The final 
metallicity averaged on the luminosity is around so- 
lar {Z « 0.01-0.04). 

The star-formation timescale is longer than, or ap- 
proximately equal to, the infall timescale of gas as- 
sembly. 

The significant majority of nearby galaxies 
(weighted by their luminosity) have Zform > 0.65, 
where Zform is the redshift that star formation be- 
gan. 

There was a peak in cosmic star-formation rate den- 
sity in the past, which was a factor of at least three 
times the present day rate. With the physical pa- 
rameterization, the peak occurs between z = 0.6 
and z = 10 (using limits of Zform < 32 and t{ > 
100 Myr). With the empirical parameterization, the 
peak occurs at z = 1 or z = 5 (with a instantaneous 
rise in the SFR at z = 5). 

Strong upper limits on star formation at high red- 
shift (a) can be obtained if we take inputs on the 



It is interesting to compute the perceived color of the cosmic spectrum to the human eye, using standard color matching functions CIE (1971, 
1986). Integrating these through the cosmic spectrum we have computed RGB values of 0.269, 0.388, 0.342. This corresponds to a blue-green 
color, the closest match in standard RGB color lists, with the same color balance (normalisation or brightness is of course arbitrary), is 'pale 
turquoise'. This is not a blackbody color, but this makes sense as the cosmic spectrum is a composite of young blue populations and old red 
populations. It is robust against choice of redshift bins. See http://www.pha. jhu.edu/'kgb/cosspec for a rendering. 
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value of low-rcdshift star formation (/3) from other 
studies (e.g. Lilly ct al. 1996). If we take /3 > 3 (i.e., 
at least 8 times the SFR at ^ = 1 relative to ^ = 0), 
then a < 0.5 (Salpeter IMF, CI). These values are 
consistent with the no dust extinction SFR models 
of Madau et al. (1998) and rule out the test case 
of a much larger SFR density at high redshift with 
a dust opacity that increases rapidly with redshift 
(their figure 7, which has /3 « 3 and a « 1). The 
test case is also not consistent with our results for 
C2 which has a < —1 for /3 > 3. This upper limit 
corresponds to a maximum of 65% and 55% of the 
stars forming at z > 1 in CI and C2, respectively. 
Correspondingly, if we take (3^2 then, a maximum 
of 80% and 75% formed at z > 1. Note that if there 
is any significant star formation prior to z = 5, this 
will lower the value of a (or /3) required to give a 
suitable model fit to the present-day cosmic spec- 
trum. 

6. Alternatively if we take the view that star formation 
was either constant or declining at high redshift (i.e. 
Q < 0) then we infer a rise in SFR to redshift unity, 
1.5 < /3 < 5, consistent with studies of luminosity 
density evolution. 

7. From the 2MASS-2dFGRS measurement of J^stars, 
Cole et al. (2001) concluded that their results were 
"only consistent with recent determinations of the 
integrated cosmic star formation if the correction 
for dust extinction is modest." Here we find again 
that there can not be massive amounts of star for- 
mation concealed at high redshift by dust extinc- 
tion. Cole et al.'s analysis compares the J and Kg 
luminosity density derived from 2dFGRS redshifts 
and 2MASS photometry with semi-analytic mod- 
els of galaxy/star-formation to derive an integrated 
stellar mass density today. Our method derives the 
star-formation history empirically using 2dFGRS 
spectra and redshifts. 

8. From the 6j luminosity density of Madgwiek et al. 
(2002) and computing the mass-to-light and SFR- 
to-light ratios from the PEGASE scenarios, we es- 
timated the present-day stellar-mass density and 
SFR density. Restricting the models to good fits 
with (3 > 1.5 and a > —3, there is still signifi- 
cant degeneracy particularly with metallicity. For 
cosmology CI with the Salpeter IMF, flgtarsh is in 
the range 0.0020-0.0062 and psFR is in the range 
0.024-0.065 h Mq yr^^ Mpc"^, and with the Kenni- 
cutt IMF, Ostarsft- is in the range 0.0013 0.0033 and 
PSFR is in the range 0.012-0.055 HMq yr-^ Mpc"^. 
For cosmology C2, the mass-to- light ratios increase 
and the values for Qgtarsh are increased by a factor 
of about 1.2. The SFR- to- light ratios remain ap- 
proximately the same. These results are consistent 
with the stellar-mass density derived by Cole et al. 
and with other derivations of the local SFR density 
(e.g. Sullivan et al. 2000). 

Finally, we compare our results with two different com- 
pilations of cosmic SFR based on the rest-frame UV 
luminosity-density technique and with a compilation of 



various other techniques. The a-/3 empirical parameter- 
ization was fitted to the data sets after converting to cos- 
mology CI. The various contours are shown in Figure 13 
with full references quoted in the caption. 

For the first UV compilation, we took the 10 data points 
with error bars from figure 9 of Steidel et al. (1999) and cal- 
culated formal confidence boundaries in a versus /?. The 
no extinction and extinction corrected data were consid- 
ered separately. This compilation promotes a high value 
oi l3 ^ 4 with a depending on extinction. 

For the second UV compilation, we took 9 data points 
with error bars from figure 13 of Cowie, Songaila, & Barger 
(1999) for the 2800A luminosity density excluding those 
used in the first compilation. A minimum error of 0.07 in 
the log was used to calculate the boundaries. This compi- 
lation promotes "a gradual decline" in the UV luminosity 
density for z < 1.5 with a = (3 ^ 1.5. 

For the other compilation, we recomputed the bootstrap 
error from the measurements of /? given by Hogg (2002) 
excluding the UV measurements included above. This new 
value was /3 = 3.2 ± 0.9. 

A precise interpretation of Figure 13 depends on the 
the accuracy of the errors used in the analysis. Given the 
analysis presented here, the two UV compilations (taking 
either extinction model from the first) are consistent at 
the 3a level with 1.8 < ^ < 2.9 and -1.0 < a < 0.7, 
but not quite consistent at the 2a level. Our results are 
independently consistent with these ranges of values but 
do not constrain the region further. Thus, evolutionary 
synthesis applied to the 2dFGRS-measured cosmic spec- 
tra {z 0.03 0.25) is in concordance with the best-fit re- 
sults obtained from a variety of rest-frame UV luminosity- 
density measurements {z ~ 0-4.5). 



6.1. Future work 

With this technique of primarily using the 'high- 
frequency' spectral information, there are still significant 
degeneracies in determining the cosmic star-formation his- 
tory. However, combined with luminosity-density meth- 
ods and with further improvements, there is the poten- 
tial to obtain more accurate star-formation scenarios and 
to discriminate between IMFs and cosmologies. For ex- 
ample, the separation of galaxies into groups with similar 
star- formation history may help break the degeneracy, i.e., 
reducing the merging of features between young and old 
populations. In addition, it may be possible to improve 
the fit to the spectra by combining models with different 
chemical evolution scenarios, e.g., weighted combinations 
of spectra with different ti.f-s- ^ibrm- 

Future large surveys such as the SDSS main galaxy sam- 
ple will provide a consistency check with the 2dFGRS spec- 
tra (selected with different effective wavelengths, 6200A 
versus 4700A, magnitude limits and aperture diameters, 
see Figure 10). With the 2dFGRS spectra (resolution 8- 
9A), and especially so with the SDSS spectra (3-4A), it is 
apparent that higher spectral resolution population syn- 
thesis models (^2A in the optical wavelength range) are 
needed to maximize the scientific potential of these data 
sets. Work is currently underway for this purpose. 

Both galaxy surveys can also be used to measure the 
slope of luminosity density variation with redshift in the 
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Fig. 13. — Comparison between different cosmic SFR studies; confidence limits in a versus f3. The solid contours represent the 2cr and 3cr 
confidence limits from this paper with POM B, r = 1.1, Salpeter IMF and Co = (0.70, 0.3, 0.7). The dashed contours represent the limits from 
a compilation of UV luminosity-density measurements with a redshift range of 0.2-4.5 (Connolly et al. 1997; Lilly et al. 1996; Madau et al. 
1996; Steidel et al. 1999), see Steidel et al.'s figure 9. The lower and upper dashed contours are for "no extinction" and "extinction corrected". 
The dotted contours show the limits from a compilation of UV measurements with a z range of 0.0-1.5 (Cowie et al. 1999; Loveday et al. 1992; 
Sawicki et al. 1997; Treyer et al. 1998), see Cowie et al.'s figure 13. The diagonally shaded region represents the la bootstrap confidence on 
/3 from the compilation of various techniques by Hogg (2002) except we have excluded the UV measurements. The data includes nebular-line 
emission, far infrared and radio continuum measurements with a z range of 0.0-1.0 (Flores et al. 1999; Gallego et al. 1995; Glazebrook et al. 
1999; Haarsma et al. 2000; Hammer et al. 1997; Hogg et al. 1998; Mobasher et al. 1999; Rowan-Robinson et al. 1997; Tresse & Maddox 1998). 



local universe which provides another star formation esti- 
mate from the same survey. This can done by either by 
using k-corrections and determining the luminosity func- 
tion in the rest frame or by determining the luminosity 
function in observed wavelength, i.e., the extra-galactic 
background light per unit redshift. These measurements 
can be used to apply further constraints on the cosmic 
star-formation scenarios and the absolute cosmic comov- 



ing star-formation rate density. 

The 2dF Galaxy Redshift Survey was made possible 
through the dedicated efforts of the staff of the Anglo- 
Australian Observatory, both in creating the 2dF instru- 
ment and in supporting the survey observations. Karl 
Glazebrook and Ivan Baldry acknowledge generous fund- 
ing from the David and Lucille Packard foundation. 
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